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This- application claims the benefit of U.S. Provisional patent application number 
60/076,698, filed March 3, 1998, the disclosure of which is incorporated herein by reference. 

FIELD OF THE INVENTION 

The present invention relates to non-invasive measurement of the mechanical 
properties of bone and especially to measurement at arbitrary angles relative to the bone. 

BACKGROUND OF THE INVENTION 

It is known in the art that the velocity of a sound wave in a material depends on the 
mechanical properties of the material. 

A sound wave that reaches a solid at an angle will typically propagate through and 
along the solid as a combination of three waves, namely, longitudinal, transverse and surface 
waves, wherein each wave has a different velocity. In bone velocity determination, the 
longitudinal wave, which is the fastest, is usually measured. The velocity of the longitudinal 
wave is: 



where E, a and p are, respectively, the Young's Modulus, the Poisson's ratio of lateral 
contraction to longitudinal extension and the mass density of the material. 

In an article entitled, "Osteoporotic Bone Fragility: Detection by Ultrasound 
Transmission Velocity," R.P. Heaney et al., JAMA . Vol. 261, No. 20, May 26, 1989, pp. 2986 
- 2990, the Young's Modulus of bone, E, is given empirically as: 

E^Kp^ (2) 
The velocity of the longitudinal sound wave in the bone is then: 

Vl-^(E/p)^4(K^) (3) 

where K is a constant which incorporates a number of factors, such as spatial orientation of the 
bone stmctures, inherent properties of the bone material and fatigue damage. Thus, the 
velocity of a longitudinal wave is a function of the mass density and can be used as an 
indicator of the quality of bone. 

In order to perform in vivo ultrasonic measurements of the mechanical properties of 
bone, it is necessary to transmit an ultrasonic wave through the soft tissue surrounding the 
bone. Unfortunately, the thickness of the soft tissue varies along the length of the bone. Also, 
the soft tissue velocity is not a constant value for all soft tissues. These variations can affect 
the accuracy of the ultrasound propagation time measurement through the bone. Typically, the 



wo 99/45348 




PCT/IL98/00299 



variations in thickness of the soft tissue and its velocity are either ignored or an attempt is 
made to cancel the effects of the soft tissue. 

U.S. Patent No. 5,143,072, the disclosure of which is incorporated herein by reference, 
describes a method of overcoming the effects of the unknown thickness of the intervening soft 
5 tissue, by ensuring that the measurements will be taken when the portion of the path which 
passes through soft tissue is of a same length for different measurements. A transmitter and 
two receivers are placed in a coUinear configuration parallel to the bone. When a wave is 
transrnitted from the transmitter towards the bone, the wave passes through intervening soft 
tissue and then travels along the bone. The two receivers detect ultrasonic waves that exit the 

10 bone and travel back through soft tissue to the two receivers. Ignoring the soft tissue, the 
difference between the path from the transmitter to the first receiver and to the second receiver 
is a segment of bone whose length is the same as the distance between the two receivers. 
Generally, the soft tissue cannot be ignored. However, if the two receivers are rather close 
together, the length of the paths in the soft tissue, between the bone and the receivers, will be 

15 approximately the same and should, to a certain level of precision, cancel out. In one 
embodiment described in the above patent document, the receiver/transmitter configuration is 
rocked and the measurements are taken only when the (shortest) distances between the bone 
and the two receivers are the same. These distances may be measured using the receivers as 
transmitter/receivers that bounce a wave off the bone. When the propagation times are equal, 

20 the configuration is assumed to be collinear with the bone. 

However, even this method has several shortcomings. First, soft tissue velocity is not a 
constant, rather, it varies with the type of soft tissue. In addition, the propagation paths 
between the bone and the receivers are not the same for the reflected wave and for the wave 
from the transmitter, so the calculated acoustic bone velocity may not be correct. Second, the 

25 above-described method requires a relatively long portion of flat bone. Thus, only a small 
number of bones can be tested, using this method, such as the tibia. In addition, since high 
frequency ultrasonic waves are very lossy, it is not practical to use them for this method. 

PCT publication WO 97/13145, the disclosure of which is incorporated herein by 
reference, describes an alternative method of bone velocity determination, in which a velocity 

30 in a significantly shorter portion of bone may be measured. In this pubhcation, several waves 
are transmitted to the bone and received by one or more receivers. One of the waves travels 
through both bone and soft tissue and one or more waves travel only through soft tissue. The 
waves that travel only through soft tissue are used to calculate the soft tissue velocity. The 
calculated soft tissue velocity is applied to extract the bone velocity from the travel time of the 

35 wave that travels through both bone and soft tissue. 
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SUMMARY OF THF. INVENTION 
It IS one of the objects of some embodiments of the present invention to provide a 
method of acoustic bone velocity detemination having a high spatial resolution. In addition, in 
some preferred embodiments of the invention a small portion of bone can be measured, so that 
ahnost all the bones of the human body can be measured. 

Another object of some preferred embodiments of the invention is to provide a method 
of acoustic bone velocity determination even when a measuring probe is not parallel to the 
bone. It should be appreciated that in many medical situations, as opposed to non-destructive 
testing situations, the relative layout of the bone and the skin cannot be accurately determined 
in a simple manner. In non-destructive testing situation, on the other hand, any intervening 
layer is usually manufactured to a known thickness. 

Another object of some preferred embodiments of the invention is to provide a method 
of acoustic bone velocity determination where the bone is not smooth and/or where the bone 
presents a curved or otherwise non-planar surface to a measurement probe. 

Another object of some preferred embodiments of the invention is to provide a 
measurement probe and/or a measurement method in which there is a minimum amount of 
mterference between transmitting elements and receiving elements. 

One aspect of some preferred embodiments of the invention relates to acquiring several 
measurements of the travel time of an ultrasonic wave, along two or more different paths 
between one or more transmitters and one or more receivers. In a preferred embodiment of the 
invention, at least three and preferably all the paths include a section of bone whose acoustic 
velocity is to be determined. 

Another aspect of some preferred embodiments of the invention relates to the amount 
of overlap between paths along which ultrasonic waves travels as a part of measuring the 
ultrasonic velocity. In a preferred embodiment of the invention, . the bone sections of the 
different paths do not substantially overiap. Additionally or alternatively, the paths of the 
ultrasonic waves in soft tissue surrounding the bone do not substantially overiap. 

Another aspect of some preferred embodiments of the invention relates to the relative 
alignments of transmitters and receivers used to deteiraine travel times for ultrasonic waves. In 
a preferred embodiment of the invention, not all of the transmitters and/or receivers are 
coUinear. In a preferred embodiment of the invention, the transmitters and/or receivers are not 
all coplanar. Alternatively or additionally, the transmitters and receivers used may be arranged 
m any order, for example, one order is a transmitter, a receiver, a transmitter and a receiver. 
Alternatively or additionally, a single ultrasonic element may function both as a transmitter 
35 and as a receiver. 
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Another aspect of some preferred embodiments of the invention relates to a method of 
solving equations v^hich links measured values and unknown variables. In a preferred 
embodiment of the invention, the acoustic bone velocity is determined by solving a set of 
simultaneous equations, into which measurement times are inserted. Additionally or 
alternatively, the simultaneous equations are analytically simplified so that only one equation, 
containing only one variable, for example the bone velocity, remains. 

Another aspect of some preferred embodiments of the invention relates to being able to 
perform several levels of estimation of a determined variable, by making different assumptions 
about other unknowns. In a preferred embodiment of the invention, at least four variables are 
interrelated: bone velocity, soft tissue velocity, an average distance of probe from bone and an 
angle between the probe and the bone. In a preferred embodiment of the invention, any one of 
these variables may be estimated or be calculated from the other variables by solving a set of 
equations. In a preferred embodiment of the invention, the soft tissue velocity is estimated, for 
example to be 1500 m/s, so only three variables need to be calculated. Thus, only three travel 
time measurements are necessary. Additionally or alternatively, more than four variables may 
be interrelated, for example, if the soft tissue is assumed to contain two layers, each with a 
different acoustic velocity. Additionally or alternatively, a known or assumed "bone" acoustic 
velocity may be used to determine a distance between a transmitter and a receiver and/or other 
dimensions of a probe, during a calibration stage with a phantom having known 
characteristics. 

Another aspect of some preferred embodiments of the invention relates to a method of 
determining travel time of an acoustic wave along a desired path. In some preferred 
embodiments of the invention, a travel time between a transmitter and a receiver are 
determined by detecting a first wave arriving at a receiver from the transmitter. In some cases 
however, that first wave may not have traveled along a desired path. In a preferred 
embodiment of the invention, a wave is determined to have traveled along a desired path based 
on characteristics of the detected wave. In one example, a spectra of a wave is different 
depending on whether it traveled through bone or not. In a preferred embodiment of the 
invention, the spectra is different due to frequency-attenuation relationships and/or frequency- 
dispersion relationships being dependent on the material through which the wave travels. 
Additionally or alternatively, if two waves arrive at a receiver, along two different paths, an 
increase in wave amplitude, at least for some frequency components, may be expected when 
the two wave arrive overlapping in time at the receiver. In a preferred embodiment of the 
invention, instead of predetermining such special characteristics of a wave, the time of arrival 
of a desired wave at a receiver is determined by performing a temporal correlation of waves at 
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two receivers which detect a similarly characterized wave. 

Another aspect of some preferred embodiments of the invention relates to reducing 
interference caused by cross talk between ultrasonic elements of a receiver/transmitter probe, 
especially for acoustic bone and/or soft tissue velocity determination. In a preferred 
embodiment of the invention, the probe design presents a labyrinth to waves which travel 
through the probe, between ultrasonic elements of the probe. Thus, such wave may be delayed. 
Additionally or alternatively, spaces between such ultrasonic elements are filled with 
acoustically slow materials and/or acoustically attenuating materials. This velocity reduction 
allows for time gating out the waves which travel through the probe. Additionally or 
alternatively, electrical shielding is provided to shield receivers from transmitters. This reduces 
electrical cross-talk. 

In a preferred embodiment of the invention, an acoustic chamber between the 
ultrasonic elements and the body to be measured is formed of attenuating materials to reduce 
cross-talk between the ultrasonic elements. Alternatively or additionally, the chamber is filled 
with acoustically slow materials. 

Alternatively or additionally, when an ultrasonic probe is placed against a human 
subject an attenuating and/or low acoustical velocity coupling material is used, to delay and/or 
attenuate undesirable cross-talk. 

Another aspect of some preferred embodiments of the invention relates to defining a 
detection window having a late cutoff time additionally or alternatively an early cutoff time. In 
a preferred embodiment of the invention, such late and/or early cutoff times are defined so that 
ultrasonic waves which do not pass through the bone cannot arrive at the receiver within the 
time window. Alternatively or additionally,, such waves may arrive, but only at an extremely 
attenuated amplitude. In a preferred embodiment of the invention, the cutoff time or times for 
each probe, are stored on a storage media associated with the probe. In a preferred embodiment 
of the invention, the storage media is packaged with the probe. Alternatively or additionally, 
the storage media is physically attached to or mounted on the probe. Additionally or 
alternatively to time windows, other calibration information and/or frequency characteristics 
information, especially for discriminating waves which travel through bone, and/or 
characteristics of the probe and/or identification information and/or usage information for the 
probe are stored on the storage media. In a preferred embodiment of the invention, the storage 
media comprises an electronic circuit embedded in the probe. Preferably, the circuit comprises 
an EPROM. In a preferred embodiment of the invention, the circuit is interrogated when the 
probe is connected into an parent device, which parent device preferably comprises circuitry 
for driving the probe and/or determining an acoustic velocity based on waves arriving from the 
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probe. 

There is thus provided in accordance with a preferred embodiment of the invention, a 
method of determining an acoustic velocity in a segment of a bone covered with a layer of soft 
tissue having an outer surface, comprising: 

determining a first travel time of a first ultrasonic wave along a first path fi-om said 
outer surface back to said outer surface which path includes at least a first part of said bone 
segment; 

determining a second travel time of a second ultrasonic wave along a second path fi-om 
said outer surface back to said outer surface which path includes at least a second part of said 
bone segment; 

determining a third travel time of a third ultrasonic wave along a third path fi-om said 
outer surface back to said outer surface which path includes at least a third part of said bone 
segment; and 

deriving said acoustic velocity in said segment of bone fi-om said three detennined 
travel times. 

Preferably, at least two of said first, second and third waves are generated 
simultaneously by a single transmitter. 

Alternatively or additionally, at least two of said first, second and third waves are 
detected simultaneously by a single receiver. Alternatively or additionally, at least two of said 
first, second and third waves each have an average fi-equency that is substantially the same, 
when generated. 

In a preferred embodiment of the invention, at least two of said first, second and third 
waves each have an average fi-equency that is substantially different, when generated. 
Alternatively or additionally, at least two of said first, second and third waves each have an 
average firequency that is substantially different, when detected. 

Alternatively or additionally, at least two of said first, second and third waves each 
have an average frequency that is substantially the same, when detected. 

In a preferred embodiment of the invention, each of said first, second and third paths 
comprises soft tissue portions and wherein at least two of said first, second and third paths 
have an overlap of at least 20% over the length of their soft tissue portions. 

Alternatively or additionally, each of said first, second and third paths comprises soft 
tissue portions and wherein at least two of said first, second and third paths have an overlap of 
at least 30% over the length of their soft tissue portions. 

Alternatively, each of said first, second and third paths comprises soft tissue portions 
and wherein no two of said first, second and third paths overlap by more than 20% of the 
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length of their soft tissue portions. 

Alternatively or additionally, each of said first, second and third paths comprises soft 
tissue portions and wherein no two of said first, second and third paths overlap by more than 
30% of the length of their soft tissue portions. 
5 In a preferred embodiment of the invention, at least two of said first, second and third 

bone parts overlap at least 20% over their length. 

In a preferred embodiment of the invention, at least two of said first, second and third 
bone parts overlap at least 40% over their length. 

- Alternatively or additionally, at least two of said first, second and third bone parts 
10 overlap at least 70% over their length. Alternatively no two of said first, second and third bone 
parts overlap by 20% or more of their length. 

In a preferred embodiment of the invention, no two of said first, second and third bone 
parts overlap by 40% or more of their length. In a preferred embodiment of the invention, no 
two of said first, second and third bone parts overlap by 70% or more of their length. 
15 In a preferred embodiment of the invention, the method comprises estimating a soft 

tissue velocity and wherein deriving said acoustic velocity comprises deriving said bone 
velocity using said estimated soft tissue velocity. Alternatively, the method comprises 
determining a fourth travel time of a fourth ultrasonic wave along a fourth path ft-om said outer 
surface back to said outer surface which path includes at least a fourth part of said bone 
20 segment and wherein deriving said acoustic velocity comprises deriving a bone velocity also 
using the fourth travel time. 

In a preferred embodiment of the invention, geometric projections of at least two of 
said acoustic wave paths onto the outer surface are parallel. Alternatively no geometric 
projections of said acoustic wave paths onto the outer surface are parallel to each other. 
25 Alternatively, said acoustic waves are generated and detected by ultrasonic elements and 
wherein said ultrasonic elements are not coplanar. 

In a preferred embodiment of the invention, said outer surface is not parallel to an outer 
surface of said bone, while said waves travel through said bone. Altematively or additionally, 
deriving comprises solving a set of simultaneous equations. 

^ ^ preferred embodiment of the invention, the method comprises repeating said 
determining of travel times and said deriving of acoustic velocity for a plurality of bone 
segments, to generate a map of acoustic bone velocity of at least a portion of a bone. 
Altematively or additionally, the method comprises repeating said determining of travel times 
and said deriving of acoustic velocity for a plurahty of orientations of travel of said waves 
35 through said bone, to generate a map of directional acoustic bone velocity of at least a portion 
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There is also provided in accordance with a preferred embodiment of the invention, a 
method of determining at least one of a set of unknowns, including an acoustic bone velocity, 
soft tissue velocity, a thickness of said soft tissue and an inclination angle of an outer surface 
of said soft tissue relative to the bone, comprising: 

determining the travel time of at least three ultrasonic waves which travel from said 
surface, to said bone, along the surface of said bone and back to said surface; and 

deriving at least one of said unknovras from said three determined travel times. 

Preferably, said at least one unknown comprises the soft tissue velocity. Alternatively, 
the method comprises assuming a value for at least one of said unknowns. Preferably, said 
assxmied unknown comprises a soft tissue velocity. 

There is also provided in accordance with a preferred embodiment of the invention, a 
probe for acoustic bone velocity measurement, comprising: 

at least four ultrasonic elements, at least one of which comprises a transmitter and at 
least one of which comprises a receiver; and 

a controller which controls said at least one transmitter to transmit at least three 
ultrasonic waves through a layer of soft tissue to a bone, which controller detects via said at 
least one receiver, at least relative travel times of said three waves, after they travel along a 
surface of said bone and which controller derives an acoustic bone velocity from said 
determined at least travel times. 

Preferably, said at least four ultrasonic elements comprise three transmitters and one 
receiver. Ahematively, said at least four ultrasonic elements comprise three receivers and one 
transmitter. Alternatively, said at least four ultrasonic elements comprise two receivers and 
two transmitters. 

In a preferred embodiment of the invention, said ultrasonic elements are not all 
coUinear. Alternatively, said ultrasonic elenients are all collinear. Alternatively or additionally, 
all of said ultrasonic elements are coplanar. Alternatively or additionally, not all of said 
ultrasonic elements are coplanar. 

Alternatively or additionally, said probe is adapted to be urged against a skin layer of a 
soft tissue and wherein said ultrasonic elements are inclined relative to said layer at an 
inclination angle. Preferably, said inclination angle is determined responsive to an expected 
acoustic bone velocity. 

In a preferred embodiment of the invention, said at least three ultrasonic waves are 
generated by a single transmitter as a single wave, which wave scatters to form said at least 
three waves. 
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There is also provided in accordance with a preferred embodiment of the invention, a 
probe for acoustic velocity determination, comprising: 

at least two ultrasonic elements, at least one of which comprises a transmitter and at 
least one of which comprises a receiver; and 

a plurality of staggered acoustic barriers which significantly attenuate ultrasonic waves 
which travel along a direct line in said probe between said transmitter and said receiver. 

Preferably, volumes in said probe between said at least one transmitter, said at least 
one receiver and said barriers comprise an ultrasonic attenuating filler. Alternatively or 
additionally, said probe comprises electrical shielding for said at least one receiver and at least 
one transmitter. Alternatively or additionally, said probe comprises at least two additional 
ultrasonic elements, which elements are also acoustically separated by said barriers. 

In a preferred embodiment of the invention, said probe is adapted for acoustic bone 
velocity determination. 

In a preferred embodiment of the invention, an acoustic chamber is defined as being 
bordered by said at least transmitter and said at least receiver and a plane to be placed against a 
body, wherein said acoustic chamber comprises acoustically attenuating material. 

There is also provided in accordance with a preferred embodiment of the invention, an 
ultrasonic probe comprising: 

at least two ultrasonic elements, at least one of which comprises a transmitter and at 
least one of which comprises a receiver; and 

an acoustic chamber defined as being bordered by said at least transmitter and said at 
least receiver and a plane to be placed against a body, 

wherein said acoustic chamber comprises acoustically attenuating material. 

In a preferred embodiment of the invention, ultrasonic elements in an ultrasonic 
velocity determination probe are comprised in an ultrasonic element grid. 

Preferably, said probe scans a bone by electrically scanning said grid. Alternatively or 
additionally, at least one of said ultrasonic elements comprises a phased array. Preferably, an 
inclination angle for said phased array elements is achieved by electrically controlling said 
phased array. 

Alternatively or additionally, said probe comprises a memory device attached to said 
probe, wherein calibration data for said probe are stored on said memory device. Preferably, 
said calibration data comprises at least one distance between ultrasonic elements. Alternatively 
or additionally, said calibration data comprises at least one vertical displacement of an 
ultrasonic element. Alternatively or additionally, said calibration data comprises at least one 
vertical displacement of a path between two of said ultrasonic elements. Alternatively or 
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additionally, said calibration data comprises a late cutoff time, after which received waves are 
ignored. Alternatively or additionally, said calibration data comprises an acoustic velocity of a 
portion of the probe between at least one of said ultrasonic elements and a surface against 
which said probe is urged in use. 

There is also provided in accordance with a preferred embodiment of the invention, a 
method of rejecting parasitic signals in an acoustic bone velocity probe, comprising: 

detecting signals arriving at a receiver, from a transmitter, ostensibly after the signal 
traveled through a portion of bone; and 

rejecting said signal if said signal arrives after a predetermined time hmit, associated 
with waves which do not travel through the bone. 

Preferably, said signal is rejected if it arrives before a second time limit. Alternatively 
or additionally, said predetermined time limit is determined based on a calibration of the 
probe. 

There is also provided in accordance with a preferred embodiment of the invention, a 
method of detecting the arrival of an ultrasonic wave from a bone, in the presence of waves 
traveling substantially only through soft tissue, comprising: 

acquiring a signal representative of said bone wave and said soft-tissue waves; and 

analyzing said signals to detect changes in amplitude in at least one frequency in said 
signal, which changes are associated with said wave from said bone. 

Preferably, said analyzing comprises determining a significant increase in amplitude of 
the signal, when a wave arrives from a bone. 

There is also provided in accordance with a preferred embodiment of the invention, a 
method of detecting the arrival of an ultrasonic wave from a bone, in the presence of waves 
traveling substantially only through soft tissue, comprising: 

acquiring a signal representative of said bone wave and said soft tissue waves at a first 
receiver; 

acquiring a signal representative of said bone wave and said soft tissue waves at a 
second receiver; and 

correlating the two signal to detect arrival of a wave from said bone. Preferably, said 
correlation is performed responsive to the generation of said wave. Preferably, said correlation 
is performed responsive to an expected arrival time window of said wave. 

There is also provided in accordance with a preferred embodiment of the invention, a 
method of acoustic bone velocity determination, comprising: 

urging an acoustic bone velocity probe onto the surface of a soft tissue layer and not 
parallel to an underlying bone surface; 
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transmitting at least one wave to said bone surface and receiving at least two waves 
from said bone surface, wherein said received waves are waves that are not reflected from said 
surface; 

measuring a travel time for each of said received waves; and 

deriving said acoustic bone velocity from said measured travel times. 

There is also provided in accordance with a preferred embodiment of the invention, a 
method of acoustic bone velocity determination, comprising: 

urging an acoustic bone velocity probe onto the surface of a soft tissue layer and not 
parallel to an underlying bone surface; 

transmitting at least one wave from said probe to said bone surface and receiving at 
least two waves from said bone surface, wherein waves are transmitted and received from 
locations in the probe, which transmission and reception locations are not coplanar; 

measuring a travel time for each of said received waves; and 

deriving said acoustic bone velocity from said measured travel times. 

There is also provided in accordance with a preferred embodiment of the invention, a 
method of acoustic bone velocity determination, comprising: 

urging an acoustic bone velocity probe onto the surface of a soft tissue layer and not 
parallel to an underlying bone surface; 

transmitting at least one wave to said bone surface and receiving at least two waves 
from said bone surface, wherein the waves connect at least three locations in the probe, each 
said location being a transmission location or a receiving location and v^herein at least one pair 
of connected locations is not collinear with any other pair of connected locations; 

measuring a travel time for said received waves; and 

deriving said acoustic bone velocity from said measurements. 

There is also provided in accordance with a preferred embodiment of the invention, a 
method of acoustic bone velocity determination, comprising: 

transmitting at least one ultrasonic wave to a bone, which wave travels along the 
surface of the bone; 

receiving said wave; and 

analyzing a travel time of a particular frequency in said received wave, wherein said 
particular frequency is related to an expected thickness of a cortex of the bone. 

Preferably, said wave is transmitted as a narrow band wave at said particular frequency. 

There is also provided in accordance with a preferred embodiment of the invention, a 
method of calibrating a probe including a plurality of ultrasonic elements including at least one 
transmitter and one receiver, comprising: 

11 
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coupling said probe to a plate having a first known acoustic velocity and measuring a 
first plurality of travel times between at least two pairs of said ultrasonic elements; 

coupling said probe to a plate having a second known acoustic velocity and measuring 
a second plurality of travel times between at least two pairs of said ultrasonic elements; and 

determining from said travel times at least two distances between pairs of said 
ultrasonic elements. 

Preferably, determining comprises determining at least one average vertical 
displacement of at least one of said pairs of elements. 

There is also provided in accordance with a preferred embodiment of the invention, a 
method of calibrating a probe comprising a plurality of ultrasonic elements including at least 
one transmitter and one receiver, the method comprising: 

coupling said probe to a plate having a furst known acoustic velocity and measuring a 
first plurality of travel times between at least two pairs of said ultrasonic elements; 

coupling said probe to a plate having a second known acoustic velocity and measuring 
a second plurality of travel times between at least two pairs of said ultrasonic elements; and 

determining firom said travel times at least one average vertical displacement of at least 
one of said pairs of elements. 

In a preferred embodiment of the invention, said plurality of ultrasonic elements 
comprises at least four ultrasonic elements. 

There is also provided in accordance with a preferred embodiment of the invention, a 
probe for acoustic bone velocity determination, comprising: 

at least two ultrasonic elements, at least one of which comprises a transmitter and at 
least one of which comprises a receiver; and 

a memory device attached to said probe, wherein calibration data for said probe are 
stored on said memory device, 

wherein said calibration data comprises at least one vertical displacement of a path 
between two of said ultrasonic elements. 

There is also provided in accordance with a preferred embodiment of the invention, a 
probe for acoustic bone velocity determination, comprising: 

at least two ultrasonic elements, at least one of which comprises a transmitter and at 
least one of which comprises a receiver; and 

a memory device attached to said probe, wherein calibration data for said probe are 
stored on said memory device, 

wherein said calibration data comprises a late cutoff time, after which received waves 
are ignored. 
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In a preferred embodiment of the invention, said calibration data comprises at least one 
distance between ultrasonic elements. Alternatively or additionally, said calibration data 
comprises at least one vertical displacement of an ultrasonic element. Alternatively or 
additionally, said calibration data comprises an acoustic velocity of a portion of the probe 
between at least one of said ultrasonic elements and a surface against which said probe is 
urged in use. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will be more clearly understood by reference to the following 
description of preferred embodiments thereof in conjunction with the figures, wherein identical 
structures, elements or parts which appear in more than one figure are labeled with the same or 
similar numeral in all the figures in which they appear, in which: 

Fig. 1 illustrates schematically various variables of an acoustic velocity determination 
measurement, in accordance with a preferred embodiment of the invention; 

Fig-^llvistrates apparatus for acoustic bone velocity determination, in accordance with 
15 a preferred embodiment of the invention; 

Fig -iis a schematic illustration of an alternative embodiment of the present invention 
utilizing an array of piezoelectric transducers; 

Fig. 4 is a schematic illustration of the array of Fig. 3, illustrating the connections of 
the transducers to control and signal processing elements; 
20 Fig. S^llustrates using a phased array acoustic probe, in accordance with a preferred 

embodiment of the invention; 

Fig. 6^ and Fig.^ illustrate an effect of ultrasonic element pitch on an operating 
characteristic of a probe in accordance with a preferred embodiment of the invention; 

Fig. 8 is a schematic cut-through view of a probe in accordance with a preferred 
25 embodiment of the invention; 

Fig._£^llustrates an ultrasonic element of the probe of Fig. 8, in accordance with a 
preferred embodiment of the invention; 

Fig. 10 illustrates a negative interaction between a non-parallel bone presentation and 
one method of ultrasonic measurement of reflections from a bone; 
^0 Fig- 1 1 illustrates an effect of non-parallel bone presentation on a measurement method 

in accordance with a preferred embodiment of the invention; 

Fig, 12 is a partial schematic view of a cut human bone; and. 

Fig. a graph showing a relationship between the thickness of an object and the 

velocity of an ultrasonic wave along its surface. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Fig. 1 illustrates a simplified situation in which an acoustic bone velocity of a bone 18 
is to be measured through an intervening soft tissue layer 22. An ultrasonic wave is generated 
by an ultrasonic element 201 (a transmitter) and detected by an ultrasonic element 202 (a 
5 receiver). There are many paths that the ultrasonic wave can and does take between the two 
ultrasonic elements. However, each path has different attributes. In particular, one path is the 
fastest path. In a preferred embodiment of the invention, the fastest path includes a bone 
segment. 

cp is defined as the angle between the line connecting the two ultrasonic elements and a 
10 line parallel to the surface of bone 18. h is defined as the average distance between the 
ultrasonic elements and bone 18. A skin 20, which defines the surface of soft tissue 22 is 
usually flush with probe 250. Considering the configuration in which (pi^O, elements 201 and 
202 are situated at posinons 205' and 209' respectively, the three main paths (for a 
longitudinal wave) are: 

15 (a) a direct path along the surface of soft tissue 22, fi-om element 201 to element 202. 

(b) a path in which a wave is reflected off bone 18 at a position 210'. 

(c) a path which travels through soft tissue fi-om position 205' to a first position 206' 
on bone 18, enters bone 18 at position 206', travels along the upper layer of bone 18 to a 
second position 207' and then exits the bone and travels through the soft tissue to position 

20 209'. 

A. Acoustic Wave propagation only via Surface of soft tissue 

This wave travels through soft tissue 22 firom transmitter position 205' to receiver 
position 209' along a surface layer of the soft fissue 22. In many cases, the wave will actually 
travel along skin 20 and/or along acoustic grease which may be used to couple the ultrasonic 
25 elements to soft tissue 22. The acoustic wave travels through soft tissue 22 at a soft tissue 
acoustic velocity Vg, which is approximately 1540 m/s. 

The acoustic propagation time of a wave traveling fi-om transmitter 201 to receiver 202 
is determined by the following expression: 

rj. _ ^2 05' -209' 

T^20y-209' — > (4) 

30 where A205'-209* is a distance between transmitter 201 and receiver 202. 

B- Acoustic Wave propagation fi-om Transmitter to Receiver by reflection fi-om Bone 

This wave travels firom transmitter position 205\ through soft tissue 22, is reflected 
fi-om the surface of bone 18 at a reflection position 210' to reach receiver 202 at receiver 
position 209*. 

14 
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The propagation time for this wave is determined using by following expression: 
^20y-2I0'-209' , (5) 

where A205*-210*-209' is the length of a path from position 205' to reflected position 210', and 
then to receiver position 209'. The distance A205'-210'-209' depends on the thickness h of soft 
5 tissue 22 and angle cp. Since the path is longer than the direct path described in "A", if the soft 
tissue velocity is the same for both paths, then T205'-210*-209' ^ T9Q5».209*- 
C. Acoustic Wave propagation from Transmitter, through Soft Tissue, along the surface of 
Bone and to Receiver 

A Critical angle y is defined by a ratio between the soft tissue acoustic velocity 
10 Vg and a bone acoustic velocity Vq: 

sm V - — 

(6) 

where Vg > Vg. 

This wave travels through soft tissue 22 from transmitter position 205' to reach first 
position 206' at the surface of bone 18; the wave then enters bone 18 at a first angle with 
15 respect to a perpendicular line to the surface of the bone 18, propagates along the surface of 
the bone 1 8 from first position 206' to second position 207*, exits bone 1 8 at a second angle 
and then travels through the soft tissue 22 until it reaches receiver position 209'. 

As is well known in optics, the fastest path between positions 205' and 209' and which 
includes the bone determines the first and second angles to be the critical angle. If no such path 
20 exists, a path of type B will be faster. It should be appreciated that even though paths of type A 
may theoretically be faster than paths of type C, if the ultrasonic wave is transmitted as a 
focused beam, substantially all the energy of the beam will be directed at the bone. 

The travel time for a wave of type C may be determined using either of two methods. 
First, the duration of the travel may be measured, under circumstances described below where 
25 it can be assured that the detected wave is a type C wave and not a type B or type A wave. 
Alternatively, the travel time may be calculated as a ftmction of other defined or derived 
variables. It is useful to breakdown the travel time into components: 

T205*-206' - the time that the wave travels (in soft tissue 22) from transmitter 
position 205' to position 206*; 
30 ^206 -207' - the time that the wave travels along the surface of bone 1 8 from position 

206' to position 207*; and 
T207'-209' - the time that the wave travels (in soft tissue 22) from position 207' to 
receiver position 209'. 

15 
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Thus, the total travel time is: 

Ttotal = T205'-206' + T206'-207' T207'-209'- (7) 
The projection of path 205'-206'-207'-209' onto a line parallel to bone 18, is indicated 
by reference A212'-213', the reference used hereafter as the length of the projection. The 
5 vertical distance from position 205' to bone 18 (at a position 21 1'), is indicated by a reference 
^205-21 r- A distance A205'-212' is half the vertical distance between positions 205' and 
209'. A205'-212' A205'.21 1 ' m^Y be determined by: 

^20y'2ir - -y sm^j? 

A " (8) 
A ^211-213' . 
'^205''2J2' = 2 ^^"^ 

The distance betv^een transmitter position 205* and first position 206* (i.e., A205'.2060 
10 may be determined by: 

, ^21Z'213' . 
h- sm^ 

^205'-206' ^ — (9) 

cos/ ^ 

Accordingly, T205*-206* be determined by: 

, ^212'-213' . 
h- sm^ 

^205'-206' = Z . (10) 

T207*-209' be determined using the following formula: 

, , ^217-213' . 
sm^ 

15 '^20r~209' = ^i7~Z (11) 

F^-.cos;^ 

T206*-207 by: 

. ^217-213 ' 2hizn(x) 
T20&-207 = ^os<p - ' (12) 

Therefore, Ttotal be calculated according by: 

_ 2h cos r , A21T'^13' 
Ttotal - y^ + cos (13) 

20 In a preferred embodiment of the invention, a plurality of measurements of T^q^i are 

performed under different conditions, so as to generate a set of equations in which Vs, Vq, (p 
and h are dependent on the plurality of measurements of Ttotal (Y is a direct fimction of Vg 
and Vb, so it is not treated as a separate variable). In some preferred embodiments of the 
invention, one or more of the above dependent variables may be estimated, thereby 

25 simplifying the equations and requiring fewer measurements of Ttotal- 
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Fig. 2 illustrates an ultrasonic measurement system 260 suitable for such multiple 
measurements, in accordance with a preferred embodiment of the invention. System 260 
includes a probe 250 including four ultrasonic elements, preferably two transmitters 251 and 
252 and two receivers 253 and 254, Preferably, the ultrasonic elements comprise piezoelectric 
5 ceramic transducers. As can be appreciated, four different transducer pairs may be defined in 
this probe, each of which provides one independent measurement of Ttotal ^r the above 
mentioned set of simultaneous equations. In other configurations, the number of transmitters 
and/or receivers may be different. Additionally or alternatively, at least one of the ultrasonic 
elements may be used both as a transmitter and as a receiver. In one exemplary embodiment, 

10 one transmitter is matched with three receivers, thereby defining only three wave paths. The 
requirement for a fourth path for solving the four equations in four unknowns is eliminated by 
estimating a value for one of the unknown variables. Alternatively, a single receiver may be 
matched with three transmitters. Additionally or alternatively, extra paths (a number greater 
than the number of variables/equations) may be defined, to provide a more robust estimate of 

15 the acoustic bone velocity. 

System 260 preferably includes a first high- voltage short-impulse driver 25 6 A coupled 
to first transmitter 251, and a second high-voltage short-impulse driver 256B coupled to 
second transmitter 252. The voltage drivers preferably translate a signal firom a controller 258 
into a high voltage signal having a peak to peak voltage of 300V. System 260 preferably 

20 includes first and second signal processors 25 7 A, 25 7B. In addition, some signal processing 
may be performed by controller 258. Each of signal processors 257A, 257B preferably 
includes an amplifier and a transforming device. Signal processor 257A is coupled to receiver 
253, and signal processor 257B is coupled to receiver 254. Controller 258 is coupled to first 
and second high-voltage short-impulse drivers 256A, 256B, and to first and second signal 

25 processors 257A, 257B. Controller 258 preferably regulates the transmission, reception and/or 
preliminary processing of acoustic waves. A computer 259 is preferably coupled to controller 
258 for fiirther processing of information received fi-om controller 258 and/or for controlling 
controller 258. Computer 259 preferably stores measurement results and/or patient information 
and/or other information in local memory, preferably a database. Additionally or alternatively, 

30 computer 259 may transmit and/or receive such infomiation fi-om a remote location, utilizing 
wired, wireless and/or computer readable media. In a preferred embodiment of the invention, 
stored results are used to compare a particular measurement to other measurements of the same 
patient, at different times and/or at different body locations and/or probe orientations. 
Additionally or alternatively, the stored results may be used as a base line for various 

35 populations. 
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In a preferred embodiment of the invention, system 260 is used in the following 
manner. As a preliminary step for a bone velocity determination process, an acoustic grease, 
preferably a water based acoustic coupling gel, silicone gel and/or mineral oil, is applied to the 
surface of skin 20. Probe 250 is then positioned on the acoustic grease. The measurement 
5 process is preferably preceded by a test measurement from which it may be determined that a 
bone is actually in the field of view of the probe and/or that a proper probe size has been 
selected. 

A plurality of measurements of wave propagation time are initiated by generating 
ultrasonic vvaves at one or both of transmitters 251 and 252. If the two transmitters are 

10 simultaneously excited, they preferably utilize different wavelengths, so that it is possible to 
detemiine, at receivers 253 and 254, which transmitter initiated an ultrasonic wave. In a 
preferred embodiment of the invention, the transmitted waves have a very fast rise time, for 
example 20 nanoseconds. Additionally or alternatively, the waves have a moderate duration, 
for example 0.5|is. In a preferred embodiment of the invention, the rise time is made as short 

15 as possible, to increase the temporal accuracy of detecting a first ultrasonic wave received 
from a bone. The duration of the wave is preferably short enough so that the probe and 
measured tissue relax between consecutive measurements and/or so that any acoustic waves 
created by the transmitters die out. 

In a prefen-ed embodiment of the invention, transmitters 251 and 252 transmit a 

20 broadband pulse, preferably centered around 1.25 MHz. It should be noted that since the 
distance that the waves travel in the bone are short, frequencies higher than those used in the 
prior art are practical, in spite of the higher attenuation of high frequency sound waves in bone. 
In general, higher frequencies give more precise results than do lower frequencies. Thus, 
higher frequencies may be used in some preferred embodiments of the invention, for example, 

25 2, 5, 10 or even 20 MHz. Typically the first wave detected by the receivers is a narrow band 
wave including mainly the frequency component which travels fastest in a particular 
configuration. Thus, a single broadband transmitted wave may be used. This wave comprises a 
plurality of frequencies, with the most appropriate frequency being automatically selected by 
the filtering properties of the soft-tissue and bone. As described below, the identification of the 

30 fastest frequency may be used to determine certain other characteristics of the bone, for 
example, a thickness of its cortex. 

In a preferred embodiment of the invention, four travel times are measured, namely: 
^25 1-253' T25i-254> T252-253 T252-254- The nomenclature Ax-y is used to indicate the 
distance between point x and point y. When x and y are ultrasonic elements, the distance is the 

35 distance between the elements. In a preferred embodiment of the invention, the travel times are 
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measured by detecting a first ultrasonic wave which arrives at an ultrasonic receiver after one 
or both transmitters are energized. 

In a preferred embodiment of the invention, a noise amplitude threshold is defined 
which a detected wave must pass. Preferably the noise threshold is a characteristic of a 
particular probe and/or of a particular receiver, transmitter and/or a receiver/transmitter pair. 
Additionally or altematively, the noise threshold may be determined by measuring noise levels 
prior to energizing the transmitters. Preferably, a time window is defined for each 
measurement. Alternative methods of determining travel time are described below. 
Additionally or alternatively, other travel time determination methods known in the ait may be 
used. The following equations may be set-up, based on the methodology discussed above, with 
respect to Fig. 1 and utilizing various known values and unknowns of system 260 and utilizing 
measured travel times: 

For a first travel time, T25 1-253 • 

^B-T 25 1.253 ' "^^^^ = ^ 251-253 -^os^ (14a) 
- /, -^251-254-^251-253 



For a second travel time, T252-254- 
^^2 

Vb • T252-254 " = ^252.25^-Cos (p (14b) 



~ h ^ "^251^254-^252-254 
A 2 - /2 + • sm (p 



For a third travel time, T251-254: 
2h 

^B • T 25 1-254 " -^^^^ = ^257-25^-cos <p , (15) 



For a forth travel time, T252-253- 

Vb • T252^253 " i;^^ = (A251^253 ^ ^252-254 ' A251-254)'<^os <p (16) 

Thus, four equations utilizing the four unknowns (Vb, Vg, h, and cp) are generated. If 
one of the unknowns is estimated, for example by estimating a soft tissue velocity to be 1540 
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m/s, only three equations are needed to solve for the three unknowns. 

In a preferred embodiment of the invention, the equations are solved for each of the 
unknowns. As can be appreciated, some solving methods may be perfomied in an off-line 
manner, so that a suitable software may be generated where the measured values are simply 
plugged into the software to discover the values of one or more unknowns. Alternatively or 
additionally, some solution methods may suggest a particular algorithmic construct. 
Additionally or ahematively, other methods of solution known in the art may be used, for 
example, iterative methods. 

Solving equation (16) for h: 
, _ ^ B'T251-254'^^(r) ^251-254- tan^r ; 

^ - - - • Qos<p (17) 

Solving for cp, using equations (15), and (16): 

cos<p = ■ ^ ^J_LA^ ^"/-^^^ (18) 

2^251^254 - ^251-253 " ^252-254 

applying equations (17) and (18) to equations (14a) and (14b), yields: 



VB-(T251-253'T251- 254) 



( A251 -254' A251 -253)' 



VB ( T252 - 253 -T251- 254) 
^ A25J - 254 ' A251 - 253 ' A252 - 254 



A251-254' A251-253 
tan 7^ 



yB('^252-253-T251-254y 



(2* A 251 -254- A251-253- A 252 - 254)^ 
Vb( T252 - 254 -T251- 254) + 



= 0 



( A251 - 254 - A252 - 254) ' 



VB ( T252 - 253 ' T251 - 254) 
^•A251- 254 -A251- 253 " A252 - 254 



A251-254- A252-254 



Vb(T252 - 253 ' T251 - 254 )' 



= 0 



(19) 



(20) 



(2 * A251 - 254 - A251 - 253 " A252 - 254 )^ 
These equations may be solved, for example by defining Vg as a function of Vg, using 
one equation to define one variable W\Xh respect to the other and then applying the result in the 
other equation yielding a single equation with only Vg as a variable, hi one example, Vg is 
determined as a function of Vg using equation (19): 
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( ^251-254-^251-253)^ ^ 

(^251^254'A251-253)^ . 


^ (^251-254-^251-253) '( ^252 - 253 '^25 J- 254 )^ 

(^251-254-^251-253)^ (2 ^251 - 254 '^251- 253 ' ^ 252 - 254)^ 



M = 



('^251-254-^251-253) + 



( ^251-254- ^251-253)(T252-253 ' T 251-254) 



n2 



^ • ^251-254 - ^251-253 ' ^252-254 

( A 251-254 - ^ 251-253 )^ - ( T 252-253 ' ^ 251-254 )^ 
(^ • ^251-254 - ^251-253 " ^252-254)^ 



(22) 



5 Then Vg can be determined from equation (20). It should be appreciated that the above 

set of equations is solved for all of the four imknowns. However, it is possible to rewrite the 
above equations there is only one equation for bone velocity which is dependent directly on 
the measured travel times. Additionally it should be appreciated that equation 21 may be 
replaced by an estimate for the soft tissue velocity. 

10 hi some preferred embodiments of the invention, one or more of the above four 

variables may be estimated using a different method than the ultrasonic measurements 
described herein. In a preferred embodiment of the invention, h is estimated using ultrasonic 
range measurement of the distance of the bone from the skin. Additionally or alternatively, h 
may be estimated using X-Ray imaging. 

15 Additionally or alternatively, cp may be estimated using X-Ray imaging. Additionally 

or alternatively, cp may be estimated by perfomiing two measurements of h, one at each side of 
probe 250. Additionally or alternatively, cp may be set to zero by rocking probe 250 and 
performing measurements when cp is determined to be zero. In a preferred embodiment of the 
invention, cp may be determined to be zero if two distance measurements on from points 

20 spaced apart on the probe determine a same propagation time to an underlying bone. 

Additionally or altematively, the soft tissue velocity may be calculated using 
triangulation methods, as described for example in WO 97/13145, the disclosure of which is 
incorporated herein by reference. In the method of this PCT appUcation, two sides of a triangle 
of paths through soft tissue are created by reflecting waves from a bone. The third side is 

25 formed by the distance between two transducers. Utilizing the three sides or two sides and a 
known angle (usually 90 degrees), assuming the soft tissue velocity is the same for all three 
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sides and knowing the travel times along two sides it is possible to express the soft tissue 
velocity as a flmction of the known travel times and the distance between two transmitters. In a 
preferred embodiment of the invention, the soft tissue velocity is determined along paths 
which overlap the paths of waves between transmitters and receivers in accordance with a 
5 preferred embodiment of the invention. It should be appreciated that the determined soft tissue 
velocity is an average of the velocities in fat and muscle. Typically layers closer to the bone 
comprise more muscle, in which the velocity may be 10% higher than in fat tissue. 

In a preferred embodiment of the invention, the above described measurements are 
repeated several times and averaged, to reduced the effect of noise. Additionally or 
10 alternatively, more than the required number of paths are measured. In a preferred embodiment 
of the invention, this results in a set of over-determined equations, the solution of which is 
known in the art. 

In the above discussion it should be appreciated that each measurement path is 
independent of the other measurement paths. Thus, there is no requirement that there be an 

15 overlap between paths. Further, there is no requirement that the transducers be coUinear or 
even coplanar. It should be appreciated that if the ultrasonic elements do not all he on a same 
plane, the effective "h" for a path defined by a pair of ultrasonic elements will not be the same 
for all pairs. Therefore, as described below, corrections for the effective "h" for each path may 
be determined. Thus, when the above equations are solved for a particular path, different *h" 

20 values are used for each path. Alternatively or additionally, a height correction constant may 
be used. Alternatively or additionally, the height correction may be determined if enough path 
travel time measurements are acquired. 

Generally however, a large degree of overlap between is desirable, to reduce errors 
caused by variations in acoustic properties of different body areas. One major source of this 

25 variability is that bone is not a homogeneous material. Rather, bone comprises several layers, 
each with a different hardness (and acoustic velocity). The thickness and composition of layers 
may vary along the length of a bone. In addition, a bone typically comprises both longitudinal 
and radial sectors, each with different hardness and acoustic velocity characteristics. In some 
cases the characteristics change abruptly and in other cases a gradual change exists. In 

30 addition, bone usually comprises anisotropy material which presents different acoustic 
properties in different orientations. The degree and/or direction of the anisotropy qualities may 
also change as a fimction of the depth of the bone. Preferably, a system in accordance with a 
preferred embodiment of the invention is used to scan a significant portion of a bone in 
different orientations and/or along significant radial and/or axial extents. 

35 To properly compare two bone acoustic velocity results derived firom two different 
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measuring sessions, the two measurements must be perfomied on the same portion of the 
bone. In particular, the location accuracy along the longitudinal axis of the bone should 
preferably be on the order of 5 millimeters in long bones, such as the tibia. This accuracy is 
easy to attain using regular positioning methods, such as marking the location with permanent 
5 marker. However, the transverse positioning accuracy should preferably be on the order of 
hundreds of microns, due to the structure of the bone. Since achieving this accuracy is 
difficult, the probe is preferably used to scan the bone for at least a certain angular range of 
transverse orientations.. When measuring the bone acoustic velocity, the probe is moved in a 
transverse direction and a plurality of bone acoustic velocities are determined. The maximum 

10 or minimum determined value is used as the reference value for comparison to bone acoustic 
velocity measurements during other sessions. Additionally or alternatively, the acoustic 
velocity of bone 18 is measured from several directions on a plane perpendicular to the bone 
axis, since the cortex of bone 18 typically has a number of different sectors, each of which has 
a different hardness and acoustic velocity. In a grid embodiment described below, a grid 

15 device is preferably utilized to detemiine a tangential positioning and/or the orientation of a 
measurement relative to the axis of the bone. Preferably such a grid device is used to 
simultaneously image the bone. 

It should be appreciated that, with some bones, such as the vertebrae, measuring softer 
sectors may be more practical than measuring harder sectors, hence the search for the 

20 minimum velocity. The minimum determined velocity typically represents travel in the softest 
sector. A minimvun determined velocity found at a later date is also in the softest sector, thus, 
the velocity measurement is repeated at the same transverse location (same sector). 

Additionally or alternatively, the acquired velocity measurements are used to build a 
transverse and/or axial velocity profile map of bone 18, which is usefiil for bone structure 

25 analysis. 

In a preferred embodiment of the invention, probe 250 is oriented to avoid straddling 
significant body structures, such as large blood vessels. In a preferred embodiment of the 
invention, large blood vessels are detected by determining Doppler shifts in the spectra of the 
ultrasonic waves which are reflected by the blood. Additionally or alternatively, such 

30 structures are detected visually. Additionally or alternatively, the ultrasonic beams are made 
significantly wider than such structures, so that their effect is reduced and/or averaged out. 

The minimal required distance of wave propagation in the bone using the above 
described methods may be as short as 5, 3, 2 or even 1 millimeters. Several resolution 
parameters should be distinguished: the precision of axial and radial placement of the segment 

35 of bone being measured; the length of the segment of bone being measured; and the precision 
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in moving the measxired area a small amount during a single session. Another important 
resolution parameter is the accuracy of the velocity measurement itself. It should be 
appreciated that using a grid type probe, as described below allows one or more of these 
resolutions to be enhanced. 
5 The preferred distance between the transmitters and the receivers depends on the soft 

tissue thickness. Using methods described herein, high resolution mapping of relatively 
uneven bones is possible. For example, such bones include the vertebra, the small bones in the 
wrist and portions of bone near joints. In addition, it is possible to measure the bone velocity 
in both longitudinal and transverse directions, since the length of the measured bone segment 

10 can be very short, even relative to a radius of a bone. In a preferred embodiment of the 
invention, the probe is made curved. Preferably, the radius of the curve is selected to conform 
to an expected radius of an underlying bone. 

The small dimensions of the minimally required bone path segment make it possible- to 
scan with a high spatial resolution, using embodiments of the present invention. For example, 

15 to measure the acoustic velocity in a portion of the cortex of a tooth, a 10 MHz ultrasonic 
pulse can be used. Due to the high frequency of the ultrasound, the probe dimensions can be in 
the order of 3 millimeters and the resolution better than 1 millimeters. 

A single measurement of a bone velocity in some preferred embodiments of the 
invention is made over a time period only 1 .5 milliseconds long, which is faster than most 

20 body rhythms. Several measurements taken along the course of a body rhythm can be used to 
measure the effect of the body rhythm on the measurement, hi some cases the body rhythm 
adds an uncertainty in determining paths and/or changes in soft tissue velocity. Alternatively 
or additionally, it is the effect of a body rhythm on Vg, Vg, cp or h which it is desirable to 
measure. In a preferred embodiment of the invention, 40 single measurements are performed 

25 and averaged over a time period of about 100ms. The average value is presented as an 
estimated bone velocity, hi a preferred embodiment of the invention, the measurements are 
analyzed to determine their distribution and/or checked for internal consistency. An average 
mapping session typically includes measuring acoustic velocity at about 100-200 points 
around the circumference of the bone. 

30 Alternatively to using a broadband transmitted wave, a single frequency pulsed wave 

may be used, since in some preferred embodiments of the invention the only aspect of the 
wave analyzed is the time of first reception of a wave. Alternatively, other, more complex 
wave forms or pulses are used and the received waves are analyzed. It should be appreciated 
that the transmitting and receiving steps of the above described process can be performed in 

35 either order and can also be performed simultaneously. Preferably, different frequencies are 
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used for each wave. Alternatively or additionally, the pulses are timed, so that no two pulses 
arrive together at a single receiver. 

A four element probe, such as described above with reference to Fig. 2, is preferably 
constructed to be less than 100 millimeters long, more preferably less than 50 millimeters long 
5 and in some preferred embodiments less than 3 millimeters long. In a specific preferred 
embodiment of the invention, the probe is 32 millimeters long. Due to the very high accuracy 
requirements from such a probe, the probe body and/or structural elements and/or filler 
materials and/or spacing elements are preferably constructed from a material which does not 
substantially expand or shrink in the temperature range of 15-40°C. Such materials typically 

10 comprise a mixture of a material which expands when heated from 15 to 40''C and a material 
which shrinks when heated from 15 to 40°C. 

In a preferred embodiment of the invention, the distance between the ultrasonic 
elements is optimized for a certain expected soft tissue depth. Thus, a typical operational 
system comprises several probes, each suitable for a different range of depths. As described 

15 below, a probe connector may include a electrical indication of its type. Alternatively, a single 
grid-type probe, as described hereinbelow, is used. 

Typically, the probe is not constructed to an absolute precision of distance between 
each ultrasonic element. Instead, a probe is constructed with a precision of approximately 0.1 
millimeter, and the exact distances between the elements are measured using a phantom. The 

20 results of the measurement, having a typical precision of better than 2fim, are stored in 
computer 259 for use in the velocity determination as described in greater detail below. Such a 
phantom preferably comprises a plastic cylinder which has a cylindrical metal core with steps 
formed along its axis embedded within the plastic. Each step corresponds to a different known 
depth of the plastic overlying the metal core. When a phantom is measured with the probe, the 

25 acoustic velocity determination equations are preferably inverted so that the distance between 
the ultrasonic elements of the probe is now a function of known depths and acoustic velocities 
and measured travel times. Alternatively or additionally to using such a phantom, a hard plate 
having a known ultrasonic velocity may be used to directly measure propagation time between 
transmitters and receivers. In a preferred embodiment of the invention, the temperature 

30 dependence of acoustic velocity in the material is also known. The thickness of the plate is 
preferably at least 20mm to be longer than the wavelength and to delay reflections from plate- 
air boundaries.. 

In a preferred embodiment of the invention, two plates having different known 
ultrasonic velocities are utilized for calibration. In each plate, at least three wave travel times 
35 are measured, between pairs of transmitters and receivers. In a preferred embodiment of the 
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invention, the travel time between the inner pair of transmitter and receiver is not measured. 
Also during ultrasonic bone velocity measurement, this inner pair is preferably not measured. 
Two plates and three measurements in each plate yield six independent equations. Six 
variables may be calibrated from these equations. The variables are preferably the three 
5 distances between the ultrasonic elements and three average heights of ultrasonic elements 
from the measurement plate, i.e., the effective thickness of an acoustic chamber separating the 
ultrasonic elements from the plate, (described below). Each average height represents an 
average height which may be used in an equation which represents traveling of waves between 
that particular transmitter/receiver pair. Alternatively to three measurements in two plates, 
10 other numbers of measurements may be performed in other numbers of plates to yield a set of 
desired equations for calibration. 

In a preferred embodiment of the invention, a probe may be mechanically configurable 
to match a particular measurement situation. Such configuration may include changing the 
distance between ultrasonic elements and/or their pitch angle. Preferably, after the ultrasonic 
15 elements are properly positioned and fixed in place, their exact locations are preferably 
detemmied by measuring the travel time in a plate as described above. 

Reference is now made to Fig. 3, Fig. 4 and Fig. 5 which illustrate aspects of a further 
embodiment useful for scanning across a section 448 of a human body, such as an arm. 

In this embodiment, a probe device formed of an array of ultrasonic 
20 transmitter/receiver cells 450 is placed onto or wrapped around section 448 or is formed into a 
sock-like element 460. The cells of array 450 are preferably formed from a piezoelectric 
material, such as a piezoelectric plastic or ceramic. Array 450 is typically acoustically coupled 
to section 448 in a standard maimer, for example, using acoustic coupling grease. 

Typically, as shown in Fig. 4, the input and output leads of each cell of array 450 are 
25 connected to an analog matrix multiplexer 451 which, in turn, is connected to a driver 452 and 
to a signal processing unit 454. Driver 452 and unit 454 are typically controlled via a 
microprocessor 455. 

Multiplexer 451 enables each cell of array 450 to be individually accessed and is 
operative to define each cell as a receiver, a transmitter, a transmitter/receiver or as non-active. 

The cells of array 450 may individually be too small to form ultrasonic transducers for 
use in prior art methods, due to the attenuation caused by long paths through bone 18. 
Therefore, a plurality of groups of cells of array 450 in desired locations were electronically 
and selectably defined to be the ultrasonic elements. In a preferred embodiment of the present 
mvention, each cell of array 450 is a separate ultrasonic element as described herein. 
Alternatively, groups of cells are defined as transducers, as shown in the prior art. However, 
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when groups are so defined, one of the preferred operation modes described below is 
preferably used. 

A first preferred method of operation is to select cells and groups of cells that 
approximate the functionality of the embodiments described hereinabove. Thus, optimal 
5 placement of transmitter/receivers can be achieved without moving ultrasonic elements and/or 



In a preferred embodiment of the invention, a two step method is used to determine the 
configuration of array 450 as transmitters and receivers. As described hereinabove, a preferred 
embodiment of the invention uses probes which are optimized for a specific soft tissue 
10 thickness between the probe and bone 18. Using array 450 to image bone 18 it is possible to 
determine the thickness of underlying soft tissue 22, before bone velocity determination. A, a 
method of optimally configuring an array of elements, according to a preferred embodiment of 
the invention, comprises: 

(a) determining the thickness of underlying soft tissue 22; and 

(^) configuring array 450 into transmitters, receivers and transmitter/receivers 
having optimal distances therebetween, which are calculated based on the determined 
thickness of soft tissue 22. 

The thickness may be estimated using a time-of-flight method which assumes a soft 
tissue velocity of 1540 m/s. Alternatively or additionally, several pre-measurements are 
20 performed to determine a desirable bone path length and/or other parameters of the 
measurement. 

A preferred method of operation maps bones and soft tissues by operating different 
cells of an-ay 450 instead of physically moving a unit comprising a plurality of ultrasonic 
units. Thus, the bone velocity at different positions and in different directions can be measured 
25 without physically moving the apparatus. 

Fig. 5 shows a multi-element probe 250', including at least one subdivided element 
287, which is subdivided into "n" sub-elements. Probe 250* preferably comprises "M" 
elements, each of which may be subdivided. In a preferred embodiment of the invention, the 
"n" sub-elements are actuated with a predetermined phase shift of a wave, to generate a beam 
30 having desired characteristics, such as direction and focusing. Using the above-described grid 
elements, it is possible to relocate a measuring zone on surface of the bone in a direction A, 
without moving probe 250'. 

In a preferred embodiment of the invention, a single grid type imaging probe is used 
for imaging, for soft tissue velocity determination and/or for bone velocity mapping. In some 
35 cases, different measurement schemes may be used for each type of measurement. 



the probe. 
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Alternatively or additionally, two different methods of measurement may be to measure a 
single property, such as bone velocity or soft tissue velocity. 

Preferably, the transmitters and/or receivers are oriented to preferentially emit and/or 
receive their waves at an angle which is the estimated critical angle. In a preferred embodiment 
of the invention, the ultrasonic waves are transmitted as pencil beams. Altematively or 
additionally, the waves are focused at the estimated or measured distance to the bone, 
preferably along a (desired) path which the waves are expected to follow. 

It is known that the velocity of ultrasound in a purely cortical bone is approximately 
3000 m/s - 4700 m/s (which corresponds to a Critical angle of approximately 19° - 31°), and 
the velocity in trabecular bones with a thin surface layer of cortical bone is approximately 
1650 mJs - 3000 m/s (which corresponds to the critical angle of approximately 31° - 69°). 

Figs. 6 and Fig. 7 illustrate an effect of ultrasonic element pitch on an operating 
characteristic of a probe in accordance with a preferred embodiment of the invention. A width 
of the first and second transmitters 251, 252, and a width of the first and second receivers 253, 
254 are preferably selected based on desired wavelength characteristics of the ultrasonic 
elements and an angle of a pitch of the transmitters 251, 252 and/or receivers 253, 254 is 
detemiined with respect to the horizontal plane. As a result of the pitch of the transmitters 251, 
252 and/or receivers 253, 254, a region 276 is defined below the pitched transmitters 251, 252 
and/or the pitched receivers 253, 254. Region 276 is also called an acoustic chamber. Region 
276 may be filled by a hquid, by an elastic material, for example, silicone rubbers, for 
example, RTV630, manufactured by General Electric Inc., US or polyurethane rubbers, for 
example EL-1 lOH, manufactured by Robner, UK or by a relatively hard material, for example, 
hard polyurethane such as U146A and a chain extender UL143 manufactured by Polymer 
Gvulot. The acoustic parameters of region 276, including one or more of a speed of sound, 
acoustic impedance and acoustic absorption, are preferably selected to optimize different 
characteristics of probe 250. 

In a preferred embodiment of the invention, probe 250 includes at least one barrier 
element 265 (shown in Figs. 6 and 7) which is preferably an acoustic attenuator. In a preferred 
embodiment of the invention, the barrier element is used to delay and/or attenuate ultrasonic 
waves which travel through the probe itself Generally, paths along which the wave travels 
faster than along the bone inclusive paths are termed "parasitic paths". 

As shown in Figs 6 and 7, the actual angle between the ultrasonic elements and the 
bone is dependent on the edge of the transmitter firom which the wave is generated. The angle 
may be a physical angle or an angle generated using a phased array. Typically, a pencil beam 
is used so that a particular inclination angle defines a range of angles yi and 72 between which 
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measurement conditions are optimal. As described above, the required range of angles depends 
on the soft tissue velocity, angle relative to the parallel of the bone and bone velocity. These 
parameters may be controlled to some extend by a user, however, many parameters depend on 
physiological characteristics which are determined by the age or sex of the patient and/or the 
5 particular body portion which is the subject of the measurement. In a preferred embodiment of 
the invention, the distances between the ultrasonic elements as used in the above equations are 
determined by the inclination angle of the elements. In general, when a low inclination angle is 
used, the upper (intemal) end of the ultrasonic element is effective and when the inclination 
angle is high, the lower (external) end of the ultrasonic element is effective. Generally, the 

10 inclination angle is selected to higher when a lower bone velocity is to be measured. For 
example, in Fig. 6 the upper end of an element and in Fig. 7 the lower end of an element. In a 
preferred embodiment of the invention, the determination of which end is used may be based 
on an estimate of the above described variables. Alternatively or additionally, an iterative 
evaluation method may be preferred, in which the results of one stage are used to estimate the 

15 values for constants in a next stage. In a preferred embodiment of the invention where a grid 
type probe is used, the piezoelectric element itself is flat, so there is no problem of upper or 
lower ends, even when the beam is inclined. 

As indicated above, one aspect of some preferred embodiments of the invention is to 
reduce parasitic waves. In a preferred embodiment of the invention, parasitic waves are 

20 reduced by placing a plurality of acoustic barriers in probe 250. Thus, any parasitic wave is 
required to travels a torturous path between transmitter and receiver elements. In a preferred 
embodiment of the invention, the barriers are placed to create a labyrinth. In addition to 
parasitic waves which travel in the probe body, some parasitic waves can travel through the 
acoustic chamber. In a preferred embodiment of the invention, barrier 265 extends through the 

25 acoustic chamber to the skin. In some preferred embodiments of the invention, barrier 265 
protrudes from probe 250, to cause an indentation in the soft tissue when probe 250 is pressed 
against the skin. 

In a preferred embodiment of the invention, probe 250 and/or the acoustic chamber are 
formed of an acoustically attenuating material. Thus, parasitic waves cannot travel through the 

30 probe with any significant intensity. By making the acoustic chamber attenuating, any 
ultrasonic wave which does not immediately exit the chamber is extremely attenuated by travel 
through it. One example of a parasitic wave which is thus attenuated is a wave which travels 
through the acoustic chamber and enters the soft tissue to "bypass" barrier 265. In a preferred 
embodiment of the invention, the acoustic chamber comprises acoustical windows so that 

35 ultrasonic waves which travel along desirable paths (at about the expected critical angle) are 
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. less attenuated than other waves. Such windows may be formed by varying the composition of 
the material of the acoustic chamber and/or by various methods of post treatment of the 
material, 

Altematively or additionally, probe 250 and/or the acoustic chamber are filled with a 
5 material having a low acoustic velocity, so that parasitic ultrasonic waves are slowed down. 

In a preferred embodiment of the invention, the acoustic chamber is filled with U-146 
Polyurethane (1920 mJs, 5.5 dB/cm), which undergoes finely bubbled phase separation when it 
polymerizes. Altematively or additionally, polyutherane EL- 11 OH (1620 m/s, 3dB/cm), by 
Robner, UK may be used as a filer material. Altematively or additionally, PWO-02 (1835m/s, 

10 15dB/cm) or PWN-01 (2010 m/s, 20 dB/cm), by CIRS USA may be used. In a preferred 
embodiment of the invention, higher attenuation materials are used when the distance between 
a transmitter and a receiver are shorter. Preferably, the path between the inner transmitter and 
receiver is not measured, so only the next shortest path needs to be taken into account. 

Fig. 8 illustrates a probe construction in accordance with a preferred embodiment of the 

15 invention. Probe 250 preferably includes acoustic barriers, as described above and/or electrical 
shielding, for isolating the piezoelectric elements from each other. In Fig. 8, each ultrasonic 
element is electrically shielded by electrical shields on either side thereof, indicated by 
reference numbers 301A-301H. Shielded cables 302A-302D are used to connect element 
drivers and other electronics to the ultrasonic elements. The shielding of the cable is preferably 

20 connected to one or both of the electrical shields surrounding each element and then to the 
bottom face (front surface) of the piezoelectric element. Altematively or additionally, an 
acoustic barrier (300A, C, E) is located between each two ultrasonic elements. Altematively or 
additionally, acoustic barriers (300B, D) are preferably located between the other acoustic 
barriers to convert the probe interior into an acoustic labyrinth. In some cases, barriers 300C 

25 and 300E may be lowered and/or extended towards the probe surface, especially if the distance 
between ultrasonic elements is shorter. In a preferred embodiment of the invention, an acoustic 
barrier is also formed surrounding all the acoustic elements of the probe and separating them 
firom an outside casing of the probe. Preferably, the electrical shielding surrounds each 
ultrasonic element firom all sides. Preferably barriers 300A-300E comprise flat plates which 

30 extend between opposite sides of the probe. The blank spaces in Fig. 8 are preferably filled 
with an attenuating material, as described above. 

Fig, 9 shows an enlarged view of an exemplary piezoelectric element used as a 
transmitter and/or a receiver, in a preferred embodiment of the invention. 

In the embodiment of Fig. 8, the propagation time of a wave generated by transmitters 

35 251 and/or 252, which travels above skin 20 and around the acoustic barriers 300A-300E, is 
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delayed because the elements in the probe 250 are designed in a labyrinth type manner. The 
dimensions of the acoustic barriers 300A-300E and of the electrical shields 301A-301H are 
illustrated below in Table I, for two preferred probe embodiments. As shown in Table 1, the 
dimensions and characteristics of transmitters 251, 252 and receivers 253, 254, and the 
distances therebetween, vary depending if the transmitter and/or the receiver are (for a 
proximal phalanx III analysis or for a radius(hand)/tibia(leg) analysis): 
Table 1 



Probe 


Al 


A2 


Angle 


F 


L 


Bl 


B2 


type 


mm 


mm 


P 


MHz 


mm 


mm 


mm 


Phalanx 


8 


2 


23 


1.25 


23.5 


3 


11.5 


radius/ 


10.5 


4 


23 


1.25 


40.5 


4 


20.3 


tibia 

















Probe type 


B3 


B4 


Dl 


D2 


HI 


H2 


H3 


H4 


H5 


H6 




mm 


mm 


mm 


rrmi 


mm 


mm 


mm 


mm 


mm 


mm 


Phalanx 


15. 


18.5 


0.5 


1.3 


15 


0.5 


12 


17 


18.5 


26.5 


radius/ 


26.1 


31-9 


0.5 


1.75 


15 


0.5 


12 


17 


19.5 


28.5 


tibia 























10 In a further embodiment according to the present invention, a body of probe 250 is 

separated from a holder (e.g., a probe holder) encasing the probe 250 with further barrier 
substantially similar composition as that of the barriers 265, 300A-300E. Thus, a wave 
generated by transmitters 251, 252 carmot travel through the probe holder to receivers 253, 254 
and arrive before a wave which travels through a bone and/or its amplitude is reduced. This 

15 further barrier preferable has a high acoustic attenuation. As such, a gas-filled (e.g., air filled) 
gap can serve as the further barrier. In addition, the further barrier can also be composed of 
sponge materials, porous materials, etc. In a prefenred embodiment of the invention, probe 250 
is left hollow, since air has a very high attenuation at these high frequencies and/or a low 
ultrasonic velocity. In a preferred embodiment of the invention, an air-filled probe may be 
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mechanically adjusted periodically, to match specific measurement situations, such as soft 
tissue depth and expected ultrasonic bone velocity. 

Table 2 shows preferred composition materials for the body of probe 250, acoustic 
chamber 276, acousric barriers 265, 300A-300E, and electrical shields 301A-30IF, and an 
exemplary hst of manufacturers (and the respective part numbers) for these elements: 
Table 2 



Element 


Composition 


Manufacture Number 


Vendor 


Body of Probe 
(250) 


P o 1 vurpf h a n p 

A v/l jr UJ. W kilCUi^ 


iji*to-ur\ ~ cnam 
extender UL 143 


Polymer Gvulot, 
Israel 


Acoustic Chamber 
(276) 


Polyurethane 


U146-0A + chain 
extender UL 143 


Polvmpr rrvnilnt 

Israel 


Acoustic Barrier 
(265, 300A-300E) 


Neoprene (i.e., 
expanded rubber) or 
Polyurethane 


CIG-3 

2529/1 4A+ isocyanate 
44V20 


Regumi 1978 Ltd., 
Israel 

Polyurethane Ltd., 
Haifa, Israel 


Electrical Shield 
(301A-301F) 


Copper Foil tape 


8271-0050-39 


Instrument 
Specialties, USA 



The devices shown in Figs. 1-9 (and the method associated therewith) are particularly 
useful when the plane of the bone 18 is not parallel to the plane of the skin 20. Fig. 10 

10 illustrates a negative interaction between a non-parallel bone presentation and one method of 
ultrasonic measurement utilizing reflections from a bone. When a reflection is determined 
from the bone, most of the energy is reflected at the incidence angle. Thus. If the incidence 
angle is zero (when the presentation is parallel), as shown for element 251, most of the energy 
is reflected back to element 251. However, if the presentation is non-parallel, as shown for 

15 element 253, most of the energy is not reflected back to the ultrasonic element. As a result, 
there is an increased susceptibiUty to noise and more difficulty in designing a receiver which is 
sensitive over several orders of magnitude. A possible end result of this problem is a difficulty 
in determining whether there is a bone surface at a certain distance from the transmitter. 
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Fig. 1 1 illustrates an effect of non-parallel bone presentation on a measurement method 
in accordance with a preferred embodiment of the invention. Instead of measuring reflection, 
some preferred embodiments of the present invention measure only the amphtude of a wave 
which enters a bone and which is emitted by the bone. The entry of the wave into the bone and 
the emission of the wave from the bone, as shown in Fig. 11, are much less sensitive to small 
changes in the angle, as might be expected if the bone surface is uneven. This is especially true 
if the entry and exit angles and not perpendicular to the bone to begin with. 

In a preferred embodiment of the invention, once a probe is calibrated the calibrated 
parameters are then stored by computer 259 into a Read Only storage device, for example a 
FLASH ROM. This Read Only storage device is preferably coupled to probe 250 and/or 
incorporated in the probe 250, so that the cahbration of a particular probe 250 is provided on 
its own Read Only storage device. When probe 250 is coupled to a measuring system using a 
connection plug, the calibration parameters of probe 250 are preferably retrieved by measuring 
system via plug terminals and are utilized to measure velocity Vg. Alternatively or 
additionally, at least some signal processing circuitry is incorporated into probe 250 itself, 
which circuitry may utilize the calibration information. Alternatively or additionally, such 
Read Only memory may include a probe identification number and/or usage information, 
which is especially useful when a measurement system is used simultaneously with more than 
one probe connected thereto. 

One or more of the following parameters are preferably calibrated and/or stored on the 
Read Only memory: 

(a) distances between elements; 

(b) amplitudes of various parasitic waves; 

(c) time limits, minimal and/or maximal, outside of which parasitic waves may arrive 
25 at a receiver; 

(d) attenuation levels of ultrasonic waves; 

(e) transmission and/or reception frequency and/or spatial response of the transmitters 
and/or receivers; 

(f) bandwidth resonance characteristics of the transmitters and/or receivers 

(g) relative thickness of acoustic chamber underlying individual ultrasonic elements; 

and 

(h) chamber material speed of sound, which may include several speeds, one for each 
operating frequency. 

Acoustic bone velocity measurement has many uses. A first use is finding fractures and 
35 Strains in bones. When a bone is over stressed or fractured (even a hairline fracture which is 
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hard to see in X-ray images), its acoustic velocity changes markedly at the locations 
surrounding the fracture. Owing to the high resolution of some embodiments of the present 
invention, fractures in the wrist bones can be identified utilizing preferred embodiments of the 
invention. 

A second use is in estimating the density of the bone and portions thereof to determine 
the loss of minerals in the bone due to diseases of the bone, osteoporosis or low-gravity 
environments. It should be noted that the velocity is dependent mainly on Young's Modulus, 
i.e., the lower the velocity, the weaker the bone. 

A third use is to chart the healing process of a broken bone. A common practice today 
is to keep the damaged bone in a cast until a predetermined period of time has elapsed. 
However, some patients require a longer or shorter healing period. X-ray images do not 
usually show enough detail to evaluate the integrity of the bone. By measuring and charting 
changes in acoustic bone velocity, a physician can more accurately estimate the state of bone 
repair. In a preferred embodiment, a small hole is drilled in the cast and the acoustic bone 
velocity is measured without removing the cast. In some patients it is advantageous to compare 
changes in acoustic bone velocities of opposing limbs. 

Another use of the invention is for measuring the thickness of the cortex of the bone. 
Fig. 12 shows bone IS having an inner core 316 and a cortex 314. The general diameter of 
bone 18 is D and the diameter of inner core 316 is d. Thus, the thickness of cortex 314 is (D- 
d)/2. 

In accordance with a further embodiment of the present invention, computer 259 
estimates the thickness of cortex 314 through utilization of an theoretically-derived and 
empirically- validated, non-dimensional curve of normalized velocity vs. normalized thickness, 
as shown in Fig. 13 to which reference is now made. A discussion of the creation of the curve 
in Fig. 13 is discussed in the book. Stress Wave.s ir .9nlid^ written by H. Kolsky, Oxford and 
Clarendon Press, 1953. Furthermore the probe 250 shown in Figs. 1-9 can advantageously be 
used for these determinations. 

The precise shape of the curve varies with the type of material being measured. 
However, it is has been determined by the present inventors that the shape of the curve is 
approximately constant for human bones. 

The velocity Vl in the curve of Fig. 13 is normalized by the velocity Vq that would be 
achieved in an infinite solid and the thickness is normalized by the wavelength, X, of the wave 
from the fransmitter. X. is, of course, determined by Vg: 

^=VB/f (23) 
where f is the frequency of the ultrasound wave. It has been determined by the inventors that 
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the curve is approximately the same whether the thickness is the thickness D (Fig. 12) of bone 
18 or it the thickness (D - d)/2 (Fig. 12) of cortex 314. The proposed explanation is that the 
when the cortex is thick relative to X, the inner portions of the bone have no effect on the 
acoustic velocity. However, when the cortex is thin relative to X, the inner portions of the bone 
affect the acoustic velocity. The inner portions of bones are usually much softer than the 
cortex, so their acoustic velocity is much lower than the cortex's acoustic velocity. Thus, if a 
higher frequency is used, velocity in a thinner cortex bone layer can be measured. In a 
preferred embodiment of the invention, the frequency used in selected to match an expected 
bone layer thickness. Alternatively or additionally, a filter is placed on the receiver. 
Alternatively or additionally, the received wave is frequency analyzed. Alternatively or 
additionally, several frequencies are tested to detennine which frequency is most suitable. 
Altematively or additionally, a broadband pulse is used and the highest suitable (which can 
travel through a thin cortex) frequency will general arrive first, with the lower frequencies 
lagging significantly behind. 

It is noted that the curve has a region 330, for relatively small velocity ratios and small 
diameter/wavelength ratios and a region 332 for diameter/wavelength ratios greater than about 
1.5 in which the ordinate is asymptotic to 1.0. 

To estimate the thickness (D - d)/2 for a bone 18, the probe makes two measurements, 
once with a high frequency wave and once with a low frequency wave. For each measurement, 
computer 259 determines the wave velocity in the bone. Altematively, in a preferred 
embodiment of the present invention probe 250 is a broadband probe and is operated only 
once. In addition, probe 250 preferably comprises frequency filters for separating received 
high frequency signals from low frequency signals. Thus, the high frequency velocity and the 
low frequency velocity are simultaneously measured utilizing the same broadband 
25 transmission. 

The response to the high frequency input wave, which has a low wavelength X, 
provides a velocity data point 334 somewhere along the region 332 from which the velocity 
Vq can be determined. The precise location of data point 334 is unknown, since the thickness 
is not yet determined. However, its precise location is generally unimportant. 

The response to the low frequency measurement provides a velocity data point 336 
somewhere within the region 330. Because the velocity Vl is known from the measurement 
and the velocity Vq may be known from the previous measurement, the location on the curve 
of the data point 336 is known. Therefore, the ratio (D - d)/(2*X) can be determined. Since X is 
known from the frequency of the transmitter and the known velocities, the thickness of cortex 
35 3 1 4 (D - d)/2 can be determiined. 
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It should be appreciated that cortex thickness measurement is more practical when 
using the present method of acoustic bone velocity determination, than when using 
conventional methods. High frequency waves attenuate rapidly when traveling through bone 
material. So, only when the path in bone 18 is short, as is possible using the present invention, 
5 are measurements using high frequency ultrasonic waves practical. Thus, in a preferred 
embodiment of the invention, the high frequency wave used is higher than in the conventional 
devices and therefore, suitable for thinner bones. 

In a preferred embodiment of the invention, probe 250 is used to determine a different 
one of the above four variables (Vg, Vg, cp and h), than Vg. In particular, it is possible to 

10 determine the velocity in soft tissue 22 using the methods and devices described herein above. 

In a preferred embodiment of the invention, a soft tissue velocity measurement may be 
made of a tissue embedded inside a second tissue. One example of such a measurement is in 
determining the acoustic velocity of a suspected cancerous lump. In a preferred embodiment of 
the invention, soft tissue velocity is determined along two paths, one which includes the lump 

15 and one which does not. If the size of the lump is known from an imager, the acoustic velocity 
in the lump may be determined by assessing its effect on the average soft tissue velocity 
measured by the methods described herein. In a preferred embodiment of the invention, soft 
tissue scanning is effected using a grid type device, as described above, so that it is easy to 
select which path a ultrasonic pulse will travel along in the soft tissue, by varying 

20 electrification of elements, rather than by moving a probe or its elements. 

In a preferred embodiment of the invention, a layer of soft tissue may be diagnosed by 
measuring the acoustic velocity, attenuation, frequency dependent attenuation and/or 
dispersion in the layer of soft tissue. Preferably, the acoustic grease used to couple the probe to 
the body is silicone oil having an acoustic velocity of approximately 1020 m/s. Alternatively 

25 or additionally, a grease is selected that does not fill skin pores. Alternatively or additionally, 
probe 250 is not pressed too hard against the soft tissue, in a manner which might compress it. 
By changing the frequency of the wave, and/or by transmitting a broadband pulse and 
determining frequency dependent responses, it is possible, according to a preferred 
embodiment of the present invention, to penetrate the soft tissue 22 with the wave, to a 

30 predetermined depth. As such, the velocity of the pulse in a specific layer in the soft tissue 22 
(e.g., skin 20, hypodermic layer, etc.) can be determined. In general, each successive layer of 
skin generally has a greater thickness and a higher ultrasonic velocity than the layer above it. 
In a preferred embodiment of the invention, skin analysis is used for dermatology, diabetes 
diagnostics and/or other endocrinology uses. In a preferred embodiment of the invention, the 

35 distance between transmitters and receivers in probe 250 is decreased so that no wave travels 
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through bone 18 faster than through the skin. Alternatively or additionally, other 
characteristics of probe 250 may be modified, including, pitch of ultrasonic elements and 
calibration settings, especially time windows for parasitic waves. These modifications are 
preferably made so that waves which travel along upper layers of the soft tissue are detected, 
while waves which travels along other paths in soft tissue and/or even in bone, are not detected 
or arrive outside a window of time in which soft-tissue waves are expected. 

As an estimate of soft tissue layer thickness through which the fastest wave will travel, 
a thickness of about the wave-length X = c/f, may be expected, where c. is the ultrasound 
velocity in the soft tissue 22 (in m/s), and f is Uie frequency (in Hz). This fastest wave 
determines, in some preferred embodiments of the invention, the effective thickness in which a 
soft-tissue measurement is made. For example, when the frequency is approximately 1 MHz 
and the velocity in the soft tissue 22 is approximately 1500 m/s, the thickness of the layer (of 
the soft tissue 22) in which the measurement is made is approximately between 1.5 and 2.0 
mm. 



It is also possible to provide additional transmitters and/or corresponding receivers to 
determine velocities in different layers of the soft tissue. As described above, the velocity in 
the soft tissue 22 can be determined by measuring the time of a fourth pulse and solving the 
equations above for the velocity in the soft tissue 22. A velocity and thickness of a second soft 
tissue layer may be determined is two additional measurements of travel time are made. Thus, 
it is possible to determine the velocity of a pulse in, e.g., a layer of fat tissue and in a layer of 
muscle tissue. 

It should be noted that many prior art methods of bone acoustic velocity determination 
use an inexact estimate for the values of soft tissue thickness and soft tissue velocity. If an 
embodiment of the present invention is used to determine more accurate values for the soft 
25 tissue thickness and velocity, these prior art methods will give more precise results. 

hi addition, measurement of soft tissue velocity is usefiil for determination of water, fat 
and muscle content of the tissue. Thus, dehydration and rehydration of a patient can be 
analyzed by measuring the soft tissue velocity, in a selected part of the patient's body, over a 
period of time. The muscle/fat ratio of the tissue can be determined if the water content of the 
30 tissue is known, or by averaging several results taken before and after the patient drinks water. 

When scanning a human female breast, the air tissue boundary can be used as a 
reflection plane. Preferably the breast is urged against a resilient form so that it does not move 
during imaging. 

hi a ftuther embodiment of die present invention, scanning is accompUshed using a cell 
35 array as described hereinabove. Preferably, the scans include scans of the same soft tissue from 
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multiple directions so that a velocity image of the tissue can be reconstructed, preferably using 
tomographic methods. 

As described above, the arrival of a wave from the bone is usually determined by that 
wave being a first wave which arrives at a particular receiver. In a preferred embodiment of the 
5 invention, a first wave firom a bone is detected even if it is not a first wave at the receiver. In 
one preferred embodiment of the invention, use is made of the fact that the wave which travels 
through bone has different characteristics than a wave which travels only through soft tissue 
and/or a wave which travels through probe 250. One example of a different characteristic 
frequency-dependent attenuation. Another example is frequency dependent dispersion. Thus, 
10 when a wave is detected at a receiver, that wave may be analyzed to determine when the wave 
changes firom a "soft tissue" type wave to a "bone" type wave. Altematively or additionally, 
when a wave starts arriving fi-om the bone, an increase, in at amplitude of least certain 
frequency components, may be detected. Since that newly arriving wave is additive to the soft- 
tissue traveling wave. 

15 In a preferred embodiment of the invention, a desired wave is detected in a method 

which bypasses and/or supplements the above considerations. This embodiment is based on 
the following observation: If a single wave is transmitted from a transmitter and is received by 
two different receivers, the received waves will have similar characteristics. In particular, the 
change in wave which occurs when a wave arrives from the bone will be similar for the two 

20 receivers, even though the background signal caused by waves which travel through soft tissue 
is different. In a preferred embodiment of the invention, the signal traces from two receivers 
are correlated to each other. It is expected that a strong correlation be detected where the bone- 
traveled waves arrive at the receivers. Determining a time delay by correlation may, in some 
cases be more accurate that first-arrived wave and may also be more robust to noise. In a 

25 preferred embodiment of the invention, the correlation is limited to a time window in which 
waves from the bone are expected. Altematively or additionally, a periodic pulsed wave is 
used, so that the effect of a correct correlation is multiplied and can be tested against a match 
to the pulsing scheme. 



30 limited by what has been particularly shown and described herein. Rather, the scope of the 
present invention is limited only by the following claims. 



It will be appreciated by persons skilled in the art that the present invention is not 
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